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SmiMARY 

The  class  of  convoluted  Pascal  distributions  obeying  certain 
moment  conditions  are  characterized  by  a system  of  differential 
equations  satisfied  by  their  probability  mass  functions.  The  result 
contains  a characterization  of  Pascal  distributions  obtained  by 
Boswell  and  Patil  (1973)  as  a special  case. 


I.  INTRODUCTION 

Boswell  and  Patil  (1973)  characterized  several  families  of 
discrete  distributions  by  differential  equations  satisfied  by  their 
probability  mass  functions.  Among  the  families  they  studied  were 
Poisson  distributions,  Pascal  distributions  and  Binomial  distributions.  * 
Samaniego  (1976)  obtained  a characterization  of  convoluted  Poisson 
distributions  which  generalizes  their  work  on  the  Poisson  distribution. 

It  is  our  purpose  here  to  extend  their  characterization  of  Pascal  dis- 
tributions to  discrete  convolutions  of  Pascal  distributions  . Characteri- 
zations in  terms  of  differentiation  with  respect  to  appropriate  parameters 
have  been  found  useful  in  the  calculation  of  expectations  (see,  for 
example.  Price  (1958)),  and  have  also  proved  useful  in  point  and  interval 
estimation  (see  Samaniego  (1976,  1977)). 


II.  THE  CHARACTERIZATION  THEOREM 

The  family  of  Pascal  distributions  is  a two-parameter  family  of 
distributions  on  the  nonnegative  integers  having  probability  mass  functions 
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f(x|n,e)  = x=o,i,...  , <i) 

where  n=l,2,***  and  0 > 1.  This  distribution  is  often  referred  to  as 
the  negative  binomial  distribution,  although  the  usual  parametrization 
is  in  terms  of  the  parameter  p = 1/9.  The  following  result  is  due  to 
Boswell  and  Patil  (1973). 

Theorem  1.  Let  {x  } be  a family  of  nonnegative  integer  valued  random 
variables  indexed  by  parameters  0 > 1,  n=l,2,***  , and  let 

0 j 0 

= ^n,l  * Furthermore,  assume  moments  of  all  orders  exist  for  every 

n and  0.  Then 

f(xln,0)  = n[f(x-lln+l,0)  - f(xlirt-l,0)]  Vn,  V0 

if  and  only  if  X(n,0)  has  a Pascal  distribution  with  probability  mass 
function  (1). 

In  generalizing  this  result  to  convoluted  Pascal  distributions,  we 
make  use  of  several  lemmas.  The  first,  due  to  Carleman,  is  proven  in 
Breiman  (1968). 

Lemma  1.  Let  be  a sequence  of  real  numbers  satisfying 

- 1^1 

lira  k *•  W 

k 

Then  there  is  at  most  one  distribution  function  F satisfying 

^k  * 


Lenma  2.  Let  Y,  Z be  Independent  nonnegative  random  variables  with 


moment  sequences  {p,,}  and  (v,3»  respectively,  satisfying  condition  (2). 


If  X = Y + Z,  then  the  moment  sequence  of  X satisfies  condition  (2). 


Proof.  Let  X,  = EX  , k^l,2,***  . We  may  write 


‘PVk-i 


i=0  ^ 


since  absolute  moments  satisfy  (see  Lukaes  (1970), 


p . 14  ) . Thus , 


1/k  . , 1/k, 


.1/k  ^ 1/k  . 1/k 

\ ^ \ + \ 


1/k 

which  implies  lim  X^'  < *. 


Lemma  3.  The  Pascal  distribution  has  a moment  sequence  satisfying 


condition  (2). 


ki 


Proof . It  is  well  known,  of  course,  that  the  Pascal  distribution  is 
determined  by  its  moments.  Condition  (2)  is  sufficient  but  not  necessairy 
for  such  uniqueness,  and  thus  we  verify  that  it  is  indeed  satisfied  here. 


5 


Since  the  Pascal  distribution  may  be  viewed  as  the  distribution  of  the 
sum  of  independent  geometric  random  variables,  it  suffices  by  T-pmma  2 
to  show  that  the  geometric  distribution  has  a moment  sequence  satisfying 
(2).  Let  the  geometric  distribution  be  parametrized  in  terms  of  p = 1/0 
so  that  its  probability  mass  function  is 

f(xlp)  = pq^  x=0,l,2,*** 


Let  ^(t)  be  the  moment  generating  function  of  a geometric  variable. 


i.e.. 


^(t)  = p(l-qe^)"^. 


One  can  easily  show  inductively  that  “^(t)  = — r ^(t)  may  be  written  as 


j=l  J 


where  ® 1 and  1 < for  j=2,...,k.  The  k*"^  moment  of  the 

geometric  distribution  is  thus  given  by 


^ J=1  J <■ 


We  thus  have 


r9‘ 


j-l  k*' 

k 4 


If  p > 1/2,  then 


If  p < 1/2, 


Thus, 


1/k 

US  < llm  S . 

“ "k-.  P P 


In  either  case,  we  have 



lim  — ; — < ® 
k 

which  is  condition  (2). 

The  k^^  factorial  moment  of  a random  variable  X will  be  denoted  by 
and  represents  the  expectation  EX(X-l)«»*(X-k+l).  The  following 
two  lemnas  are  easily  proven  inductively. 
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Lemma  4.  Let  Y,  Z be  independent  random  variables  and  let  X = Y + Z, 

th 

Let  l^(ic)  '^(k)  ^ factorial  moments  of  X,  Y and  Z 

respectively.  Then 


’^(k)  = ^fQ(P^^(k-i)^(i)‘ 


Lemma  5.  If  Y has  a Pascal  distribution  with  parameters  n and  0, 
the  k*"^  factorial  moment  of  Y is  given  by 


(nfk- 1)1 

^(k)  (n-l)i  * 


ow  in  a position  to  prove  the  following  characterization 


Theorem  2.  Let  fx  be  a family  of  nonnegative  integer  valued  random 

variables  indexed  by  positive  integers  n and  by  0 > 1,  and  let  Z be 

a nonnegative  integer  valued  random  variable  whose  distribution  does  not 

depend  on  0.  Suppose  the  moment  sequence  <3f  Z satisfies  condition  (2), 

k 

and  EX  -♦  v,  as  0 -*  1 for  all  n,  k.  The  probability  mass  function  of 

cC 

X . satisfies 


— f(x|n,0)  = n[f (x-l|nfl,0)  - f(x|n+l,0)]  Vn,  0 


if  and  only  if  the  distribution  of  X » is  the  convolution  of  the  Pascal 

II9  w 

distribution  with  parameters  n,  0 and  the  distribution  of  Z. 


Proof . It  is  easy  to  verify  that  the  probability  mass  function  of  a 

convoluted  Pascal  distribution  satisfies  (3).  We  proceed  immediately 

to  the  converse.  Let  E(X  represent  the  factorial  moment  of 

X It  is  clear  that  since  factorial  moments  are  linear  functions 

k k 

of  the  central  moments,  EX  -*  EZ  as  9 -•  1 implies  that 

ri)0 

E(X  £,)/•.  N “•  EZ,,  , as  0 -*  1.  Moreover,  if  a probability  distribution 
n,w 

is  determined  by  its  moments,  it  is  determined  by  its  factorial  moments. 


^(k)  ^ x(x-l)*»*(x-k+l) 


and  note  that 


■ ^(k)  " ^(k-1)* 


If  {x  -}  is  a family  of  random  variables  whose  probability  mass  functions 
Q ) o 

satisfy  the  system  of  differential  equations  (3),  and  if  moments  of  all 
orders  exist,  as  postulated,  then  from  (4)  we  have  for  k > 1 


00  00 
= n|^  E (xfl)^j^jf(xln+l,0)  - E x^j^^f (xln+-l,6)] 


We  show  inductively  that  the  factorial  moments  of  X^  ^ are  those  of  a 


convoluted  Pascal  distribution  which  is  determined  by  its  moments.  Consider, 
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for  fixed  n,  the  first  moment  of  X We  have  from  (5) 

n,0 


35  - '>• 


Equating  indefinite  integrals  of  terms  in  (6)  yields 


E(X^^q)  = n(e-l)  + c. 


Since  lim  E(X  ) -*  EZ,  we  have 


0 ->  1 


E(X^^q)  = n(0-l)  + EZ. 


Thus,  for  each  n,  the  first  (factorial)  moment  of  X ^ is  that  of 

Y + Z where  Y has  a Pascal  distribution  with  parameters  n,  0. 

XI 

Assume  that  for  every  n,  the  factorial  moment  is  that  of 

Y - + Z,  that^is,  assume,  using  Lemmas  4 and  5,  that 


Consider  the  (k+l)-st  factorial  moment.  Equation  (5)  implies  that 


^ il(k-i)l  nl  i+1  ^(k-i)  ^ 


bution  and  the  distribution  of  Z. 

Remark  1 . The  characterization  theorem  for  convoluted  Pascal  distributions 
reduces  to  the  result  obtained  by  Boswell  and  Patil  when  the  random 
variable  Z is  degenerate. 


Remark  2.  It  is  interesting  that  the  system  of  differential  equations  (3) 
are  satisfied  by  any  convoluted  Pascal  distribution,  irrespective  of  the 
existence  of  moments.  The  moment  conditions  in  the  hypothesis  of  Theorem  2 
are  needed  in  our  proof  for  unique  identification  of  convoluted  Pascal 
distributions  by  their  moment  sequences.  We  conjecture  that  the  system 
(3)  alone  characterizes  the  class  of  all  convoluted  Pascal  distributions. 

Remark  3.  Let  Xj^,...,X^  be  a random  sample  from  a convoluted  Pascal  dis- 
tribution. Let  the  likelihood  function  be  denoted  by 


L(Xi , . . . ,x.  ,0)  = n f(x. |n,0) 
i=l  ^ 


. 
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